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CNJ ^ Abstract 

The generalized Chaplygin gas model (GCGM) contains 5 free parameters that must be constrained 
^ . using the different observational data. These parameters are: the Hubble constant Ho, the parameter 

' A related to the sound velocity, the equation of state parameter a, the curvature parameter ^Iko, and 

the Chaplygin gas density parameter Q^o- The pressureless matter parameter ^ImO may be obtained as a 
1/^ ■ dependent quantity. Here, these parameters are constrained through the type la supernovae data. The 

"gold sample" of 157 supernovae data is used. Negative and large positive values for a are taken into 
, account. The analysis is made by employing the Bayesian statistics and the prediction for each parameter 

is obtained by marginalizing on the remained ones. This procedure leads to the following predictions: 
5 ; " = -0-75tt% Ho = 65.0011:^^, n^o = -0.77t^it, ^^mo = 0.00lJ:gg, Qao = l-SGtHl, ^ = i.ooolg:°°°. 

, Through the same analysis the specific case of the ordinary Chaplygin gas model (CGM), for which 

^ l' a = 1, is studied. In this case, there are now four free parameters and the predictions for them are: 



Ho = 65.0ltl:n, ^ko = -2.73tlir, = 0-OOtlfo, ^co = 1.34+^:?o, A = 1.000+^;™^. To complete 
the analysis the ACDM, with its three free parameters, is considered. For all these models, particular 
^ cases are considered where one or two parameters are fixed. The age of the Universe, the deceleration 

, parameter and the moment the Universe begins to accelerate are also evaluated. The quartessence 

scenario, that unifies the description for dark matter and dark energy, is favoured. A closed (and in some 
cases a flat) and accelerating Universe is also preferred. The CGM case q = 1 is far from been ruled out, 
, and it is even preferred in some particular cases. In most of the cases the ACDM is disfavoured with 

jjj ■ respect to GCGM and CGM. 

1 Introduction 



One of the most intriguing observational recent results in cosmology concerns the possible accelerating phase 
of the Universe today. This result comes from the observational programs of the type la supernovae, carried 
out since the second half of the last decade ^E]. Type la supernovae seem to be excellent standard candles: 
their detonation mechanism occurs under very specific conditions and their absolute magnitude can be, in 
principle, easily inferred. For a more detailed discussion on this problem, see reference 3 . The crossing of 
these results with those coming from the anisotropy of the cosmic microwave background radiation ^ leads 
to a scenario where the matter content of the Universe is composed essentially of 70% of an unclustered 
component of negative pressure, the dark energy, and 30% of a clustered component of zero pressure, the 
cold dark matter. The most natural candidate for dark energy seems to be the cosmological constant, since 
it can be connected with the vacuum energy in quantum field theory ^ . But, the small value resulting 
from observations for the energy density of the cosmological constant term, pA = 10~*^GeT^*, leads to 
a discrepancy of about 120 orders of magnitude with the theoretically predicted value [B]. Hence, other 
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possibilities have been exploited in the literature, like quintessence, a kind of an inflaton self-interacting field 
adapted to the present phase of the Universe [3 El- More recently, the Chaplygin gas model (CGM) has been 
evoked [HI ■ This model is based on a string inspired configuration that leads to a specific equation of state 
where pressure is negative and varies with the inverse of the density This model has been generalized, 
giving birth to the generalized Chaplygin gas model (GCGM), where now the pressure varies with a power 
of the inverse of the density JI]. These proposals have many advantages, among which we can quote the 
following: in spite of presenting a negative pressure, the sound velocity is positive, what assures stability 
|12| : these models can unify the description of dark energy and dark matter, since the fluid can clusters 
at small scale, remaining a smooth component at large scales the CO has an interesting connection 
with string theory ^U). Some criticisms have been addressed to the GCGM (CGM) mainly connected with 
its features related to the power spectrum for the agglomerated matter ^^l- However, in our opinion, this 
specific criticism is not conclusive, since the introduction of baryons may alleviate the objections presented 
against the cosmological scenarios based on the GCGM (GCM) [Tl] . 

In order to test these different models for dark energy (and dark matter) it is essential to compare 
their predictions against the different observational results available now. These observational results refer 
mainly to the type la supernovae, the anisotropy of the CMBR, the X-ray data of clusters and super-clusters 
of galaxies, the gravitational lenses and the matter power spectrum. Each observational test constrains 
differently the free parameters of the models; the crossing of the informations coming from those observational 
data can strongly restrict or support a specific cosmological model. In carrying out this task, a statistical 
analysis must be applied for each particular observational test. One important point is how to perform this 
statistical analysis: the final conclusions may, in some cases, depend on the statistical framework (Bayesian, 
frequentist, etc.), as well as on the parameters that are allowed to be free, and how these parameters 
are constrained (through a joint probability for two parameters, minimizing the error function or through 
a marginalization of all parameters excepted one, etc.). In some cases, the different procedures adopted 
may lead to quite different conclusions on the best value for a given set of parameters. The choice of the 
observational data sample may of course be important as well. 

In preceding works |15[ [TC] we have tested the ordinary CGM and the GCGM against the type la 
supernovae observational data. The comparison between theoretical prediction and observational data, in 
what concerns type la supernovae, is made essentially by computing the luminosity function Dl and re- 
expressing it as the distance moduli 



in terms of which the observational data are given. In those works, we employed a Bayesian statistical analysis 
considering all possible free parameters, which were in number of 4 for the CGM and 5 for the GCGM. In 
reference ^SIj the same analysis has been also performed for the model with cosmological constant and 
cold dark matter (ACDM), in order to allow a proper comparison between the different models. The final 
estimation for each parameter was obtained through the marginalization on all other parameters. One 
limitation of these previous works was the use of a sub-sample of 26 supernovae data. This sample seems 
to be very restrictive when we remember that there is today the "gold" sample, with 157 supernovae, and 
the "silver" sample with 192 supernovae ^7]. Samples of up to 230 supernovae are now available But 
the 26 supernovae of those works have a very good quality, and they lead to a quite small value for the 
fitting parameter: using the "gold" sample the value for increases a little. This reflects, for example, the 
fact that those larger samples contain supernovae whose observational status is not very well established: 
for example, they contain supernovae with the almost the same redshift z but with different values for the 
luminosity distances, without a superposition of the error bars. However, the recent works on the subject 
employ enlarged samples of supernovae data for obvious reasons: they contain supernovae with 2; > 1, leading 
to a better estimation of the deceleration parameter today go; moreover, it is expected that the dispersion 
in the parameters estimations can be narrowed when more supernovae are used. Hence, in this work, we 
return back to the problem of estimating parameters but now using the "gold" sample. We will do it for 
the GCGM, the CGM and for the ACDM model. The statistical method will be the same as in references 



The type la supernovae data have been used in many works to constraint the parameters of the GCGM 
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QSj^PO]- In these works, however, restrictions on some free parameters were introduced by, for example, 
fixing the curvature of the spatial section equal to zero or using a specific value for the Hubble constant 
as suggested by the spectrum of the anisotropy of the cosmic microwave background radiation (CMBR). 
Moreover, a more simplified statistical treatment has been frequently employed, like the statistics. In 
almost all of them the parameter a has been restricted to the interval < a < 1. This last restriction has 
been alleviated in a recent works EZl EH EH HOI • For example ref. EHj , where the statistics and 
quartessence scenario (that unifies dark matter and dark energy into a single fluid, the Chaplygin gas, so the 
dark matter content is null) have been employed, values of a greater than 1 have been allowed. The best 
value found for the parameter a is 3.75 for the flat case and 2.87 for the non-flat case. 

The main reason to treat again this problem is to give a complete analysis of the GCGM, as well as of 
the CGM and ACDM, using all free parameters. The gold sample with 157 supernovae will be used. The 
parameter a will be allowed to take any positive or negative value. The Bayesian statistics will be employed, 
and the predicted values for each parameter will be obtained through the marginalization process. This last 
step leads to results quite different from those obtained through, for example, the statistics. In fact, even 
if we agree with the authors of references [211 EOI that the minimization method leads to a positive (in our 
case we find a = 5.66 when all five free parameters are free), the marginalization on the other parameters 
indicates that the best value for this parameter is negative, although the dispersion is quite large and high 
positive values can not be excluded. 

In treating three theoretical models (GCGM, CGM and ACDM) we intend to give a more unified de- 
scription of all these different cases. In what concerns the GCGM, the free parameters to be estimated 
are the value of the Hubble constant today Hq, the parameters a, A (connected with the sound velocity 
of the fluid), the curvature parameter and the GCG parameter f2cO (alternatively, the ordinary matter 
parameter rimo)- In the CGM, the parameter a is fixed to unity. For the ACDM model the parameters are 
Hq, ilco (i-e., ^Ao) and (or the curvature parameter ilto)- In both cases, we also estimate the age of 
the Universe the deceleration parameter qq and the value of the scale factor at the moment the Universe 
begins to accelerate at (which can also be expressed in terms of the redshift Zi). 

We will exemplify the subtleties connected with the statistical analysis by displaying the analysis based 
only on the better value, on the joint probability of two different parameters and by marginalizing on 
all parameters excepting one. In order also to better compare with the literature, we will also consider 
some particular cases where the spatial curvature is fixed as flat, or the ordinary matter is fixed equal to 
^niQ = 0.04 (inspired on the nucleosynthesis results) or flmo = or fixing two of the parameters. 

Our results indicate the following values for the parameters for each model. For the GCGM we find: 
a = -0.75t^;04^ Ho = 65.00^};^, fifco = -0-77tlli n^o = O-OOtUl n,o = IMtHl ^ = 1.000tl]°;^°. 
For the CGM (a = 1) the estimations give: Hq ^ 65.011^;^^, ^ko = -2.731^;^^, = O.OOtJ;^^, O^o = 
1.34+°:^;5, A = 1.00+°:°?o. Finally, for the ACDM the results give: Hq = 65.00+^;^^, Qko = -2.12^^;^^, 
^mo — l-Ollg g5, ^co = l-36to'78- Hencc, our analysis indicate that the traditional Chaplygin gas model 
can not be ruled out, at least in what concerns type la supernovae data, and that the ACDM case is not 
the preferred one. However, the dispersions are large enough so that no definitive conclusion can be made. 
Only through the crossing with other tests a more restrictive scenario can come up. However, in crossing 
the different tests, a uniform statistical procedure must be used. This work intends to be the initial step in 
this program. 

This paper is organized as follows. In next section, the model and the different relevant quantities are 
set up. In section 3 we present the parameter estimation. In section 4 we present our conclusions. 



2 Definition of the models and of the relevant quantities 

The GCGM is obtained through the introduction of a perfect fluid with an equation of state given by 

P=-4 ' (2) 
P 

where A and a are constants. When a = 1 we re-obtain the equation of state for the CGM. Henceforth, we 
will use mainly the term generalized Chaplygin gas model, keeping in mind that when a ~ 1 we have the 
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traditional Chaplygin gas model. In principle, the parameter a is restricted in such a way that < a < 1. 
However, we will allow a to take negative and large positive values. Negative values and values greater than 1 
can be potentially dangerous since they can lead to an imaginary sound velocity and a sound velocity greater 
than the velocity of light, respectively. But, these problems appear more dramatically at perturbative level. 
In this case, however, a fundamental description for the GCG must be employed (for example, using self- 
interacting scalar fields) what avoids those drawbacks. The supernovae data tests mainly the background 
model, and in this sense to enlarge the possible values of a does not bring any difficulty. 

In the GCGM we introduce also pressureless matter in order to take into account the presence of baryons 
in the Universe and also in order to verify if the unified scenario (where no dark matter is present) is favoured 
by the data. Hence, the dynamics of the Universe is driven by the Friedmann's equation 

aV k _ SnG 
a I 3 



Pm + Pc] , (3) 



a 

Pm + ^-Pm = , (4) 

Pc + 3^(p.-^) = , (5) 

where and pc stand for the pressureless matter and Chaplygin gas component, respectively. As usual, 
A; = 0, 1, —1 indicates a flat, closed and open spatial section. 

The equations expressing the conservation law for each fluid H4I5|I lead to 



1/(1+") 

j3 ' '"'^ ^1 ^3(1+0) r ' (6) 



PmO ] A , ^ 



The value of the scale factor today is taken equal to unity, oq = 1. Hence, pmo and pco =s A + B 



1/(1+") 

are the pressureless matter and GCG densities today. Eliminating from the last relation the parameter B, 
the GCG density at any time can be re-expressed as 



I- A 



1/(1+") 



Pc = Pco^A+-^^j , (7) 

where A = A/pco- This parameter A is connected with the sound velocity for the Chaplygin gas today by 
the relation -y^ = aA. 

The luminosity distance is given by [21E21 

dL = ^n , (8) 

a 

ri being the co-moving coordinate of the source. Using the expression for the propagation of light 

and the Friedmann's equation (|3Jl , we can re-express the luminosity distance as 

dL^il + z)S[f{z)] , (10) 

where 

S{x)^x (fc = 0) , S'(a;)=sinx (fc = 1) , ^(a;) = sinhx (fc = -1) , (11) 



and 



■^^''^Ifo Jo {^rrAz' + 1)3 + ^icoi^ + + 1)3(1+")(1 - A)]l/(1+") - n,o{z' + m'/^ ' 
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with the definitions 

_87rGp„io „ _87rGpco n _ ^ nQ^ 

such that the condition f2mO + f^co + f^fco = 1 holds. The final equations have been also expressed in terms 
of the redshift z — —1 + ^. 

The age of the Universe and the value of the decelerated parameter 90 = ^ ff are given by 



2 



(14) 



to 

T 7o (1 + Z'){^MZ' + If + ^co[A + {Z' + 1)3(1+«)(1 - A)]l/(l+a) - f7fco(2' + 1)2}!/ 

a«o + f^co(i - 3A) 

where T = (100/iJo) x 10^°, so that to has units of years. 

The value of the scale factor signs the start of the recent accelerating phase of the Universe, it is given 
by the roots of the equation 

-d = a{H + H^) , (16) 

and is related to the redshift value Zi such that — = tt— or = ~1 + 

Following the same lines sketched above, we can obtain the corresponding expressions for the ACDM 
model. However, it is easier just to insert in the above relations the condition ^ = 1 in order to recover the 
ACDM case. 

In order to compare the theoretical results with the observational data, we must compute the distance 
moduli, as given by relation Q). A crucial aspect of the present work is the employment of the Bayesian 
statistics, which will be outlined in the next section. The first step in this sense is to compute the quality of 
the fitting through the least squared fitting quantity defined by 

i ""mo, J 

In this expression, /ig j is the measured value, (Iq j is the value calculated through the model described above, 
'^fio i measurement error and includes the dispersion in the distance modulus due to the dispersion in 

galaxy redshift due to peculiar velocities, following ref. [T7| . 



3 Parameter estimations 

To constraint the five independent parameters for the GCGM, the four independent parameters for the 
CGM and the three parameters for ACDM, we use the Bayesian statistical analysis. The method and its 
motivation are described in detail in ref. J 51- Since there is no prior constraint, the probability of the set of 
distance moduli ^0 conditional on the values of a set of parameters {ai} is given by a product of Gaussian 
functions: 

{^^0,^ - 



p{M{o-i}) = n / ^ =exp 



(18) 



This probability distribution must be normalized. Evidently, when, for a set of values of the parameters, the 
is minimum the probability is maximum. This is a valuable information but is not enough to constraint 
the parameters. In table 1 the values of the parameters for the maximum probability (minimum x^) is given, 
using the gold sample of supernovae, for the GCGM with five free parameters and for other cases where 
the baryonic, the curvature or both are fixed. The same estimations are presented in tables 2 and 3 for the 
CGM and ACDM. Note first that the minimum values for using the "gold sample" are higher than the 
corresponding ones using the restricted sample of 26 supernovae ^1^]. Note also that from this point of 
view, the best value for the parameter a is much bigger than 1. These results must be compared with a 
more complete analysis to be presented below. For more details about the minimization process, see refs. 

[13 un. 
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Table 1: The best- fitting parameters, i.e., when is minimum, for each type of spatial section and matter 
content of the generalized Chaplygin gas model. Ho is given in km/ Mpc.s, A in units of c, to in Gy and at 
in units of oq. 





CGM 


CGM : 
fc = 


CGM : 

^mO = 


CGM : 
n^o = 0.04 


CGM : /c = 0, 

^mO — 


CGM : /c = 0, 
n,no = 0.04 




1.1107 


1.1141 


1.1107 


1.1108 


1.1141 


1.1147 


Ho 


65.08 


64.73 


65.06 


65.06 


64.73 


64.70 




-0.38 





-0.36 


-0.38 










0.00 


0.00 





0.04 





0.04 




1.38 


1.00 


1.36 


1.34 


1 


0.96 


A 


0.761 


0.811 


0.764 


0.778 


0.811 


0.834 


to 


13.29 


14.01 


13.34 


13.32 


14.01 


14.03 


qo 


-0.888 


-0.717 


-0.876 


-0.872 


-0.717 


-0.701 


0.1 


0.734 


0.699 


0.732 


0.732 


0.699 


0.695 



Table 2: The best-fitting parameters, i.e., when is minimum, for each type of spatial section and matter 
content of the traditional Chaplygin gas model. Ho is given in km/Mpc.s, A in units of c, to in Gy and at 
in units of ao- 

From the probability distribution H18|l . a joint probability distribution for any subset of parameters 
can be obtained by integrating (marginalizing) on the remaining parameters, see refs. [151 116j . This is a 
valuable information. But, in order to estimate properly a single parameter, the probability distribution 
must be marginalized on all other parameters. This in general gives a quite different result if we try 
to estimate the parameter in a two or three-dimensional parameter space. The reason is that, in such 
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Table 3: The best-fitting parameters, i.e., when minimum, for each type of spatial section and matter 
content of the ACDM model. Ho is given in km/Mpc.s, to in Gy and in units of oq. 

multidimensional parameter space, if a parameter has a large probability density but in a narrow region, 
the total contribution of this region may be quite small compared to other large regions which have small 
probability: in the marginalization process, this kind of high PDF region contributes little to the estimation 
of a given parameter. These features are exemplified in figure Q where the probability distribution function 
(PDF) for the set of parameters {a,Ho,A) in the GCGM case is displayed, with a clearly non-Gaussian 
behaviour. Hence, in what follows the estimation of a given parameter will be made by marginalizing on all 
other ones. 

3.1 Estimation of a 

In the case of five free parameters, the procedure described above gives a = — 0.75lg'24. Note that this 
prediction differs substantially from that extracted from the minimization of x^, which gives a large positive 
best value a, instead of a negative a when the marginalization is made. Concerning the best value when the 
marginalization is made, even if the best value is negative, the dispersion is quite high, so even large positive 
values are not excluded, at least at 2a level. Comparing our results with other ones already published requires 
some care due to the fact that usually in the works already quoted some parameters are fixed. Moreover, 
the allowed values for a are generally restricted to the interval [0, 1], except in references [191 1271 1^ I29L 130) . 
In reference the authors considered — 1 < a < 1, a flat spatial section and they have fixed the value of 
the Hubble constant. They found a best value around a ~ 0.4. 

If we compare with the results of ref. , the use of a large sample of supernovae has slightly diminished 
the best value, but at same time, it has reduced considerably the dispersion. Imposing that the space is flat, 
or fixing the pressureless matter equal to zero or 0.04, lead to a positive best value for a; the dispersion is also 
affected. For example, with the imposition that k = and fim = 0.04, the estimation leads to a = 1.52^^ gy. 
In figure 01 the PDF for a is displayed for the case of five parameters, and also for three particular cases, 
where one or two parameter is fixed. Note that as the baryonic matter or the curvature (or both) is fixed, the 
maximum of the PDF is displaced to positive values. In general grounds, we can state that even relatively 
large values for a are not excluded. Of course, for the CGM the value of a is fixed to unity. However, from 
the analysis for the GCGM it can be inferred that the probability to have a = 1 is 19.58% when the five 
parameters are free. Restricting to null curvature or fixing the pressureless density increases considerably 
this value: for example, when the curvature is null the probability to have a = 1 is about 90%. At the same 
time, the probability to have a > 1 is 41.98% when all parameters are free, but this value can double when 
one or two parameters are fixed. 
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Fi gure li The graphics of the joint PDF as function of (a, Hq^ A) for the generaUzed Chaplygin gas model when fixing k = 
and QmO = 0.04. It is a 3D density plot where the function (here the PDF) is rendered as semi-transparent colourful gas, 
or we can say that the plot is made by voxels (volume elements) analogous to pixels (picture elements). Where the PDF is 
minimum, the transparency is total and the colour is red, where the PDF is maximum then the voxel is opaque and the colour 
is violet; mid-range values are semi-transparent and coloured between red and violet (red, orange, yellow, green, blue, violet). 
The parameter ranges are —1 < a ^ 9, 62 ^ Hq ^ 68 and 0.5 ^ A ^ 1, with resolutions 0.25 X 0.5 X 0.01, respectively. The 
joint PDF of [a, Ho, A) for other cases (5 free parameters, 4 free parameters, etc.) show a similar 3D shape. 

3.2 Estimation of A 

In general the results indicate that the value of A is close to unity. When five free parameters are considered, 
the marginalization of the remaining four other parameters leads to A — I.OOOIQ534. This may suggest the 
conclusion that ACDM (A = 1) model is favoured. However, the accuracy of the computation, due to the 
step (between 0.01 and 0.02) used in the evaluation of the parameter, does not allow this conclusion. Instead, 
it means the peak happens for 0.98 < A ^ 1. 

In fact, fixing the curvature or the pressureless matter, the preferred value differs slightly from unity: 
with fl„i = 0.04 and fc = 0, for example, we have A = 0.987lg;^72. The situation is essentially the same for 
the CGM: when all other parameters are marginalized, the best value is unity, but with a small dispersion, 
A — 1.000;to'27o; however, this best value becomes smaller when one or two parameters are fixed. This 
occurs in both for the GCGM and CGM. In figures H and H the PDF for A is displayed, both for the GCGM 
and the CGM, in the case where the marginalization is made in all other parameters, and when the matter 
density and/or the curvature is fixed. Fixing one or two parameters displaces the maximum of probability 
in the direction of smaller value of A. This effect is more sensible for the CGM: in general, the best value 
of the CGM is smaller than in the GCGM, and the dispersion is also smaller. Note that the probability to 
have A 1 is zero when all parameters are free, both for the GCGM and CGM. But, this probability can 
become as high as 100% when one or two parameters are fixed. 

In figure|21the joint probabilities for a and A are displayed, with a non-Gaussian shape. This figure is an 
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Table 4: The estimated parameters for the generalized Chaplygin gas model (GCGM) and some specific 
cases of spatial section and matter content. We use the Bayesian analysis to obtain the peak of the one- 
dimensional marginal probability and the 2cr credible region for each parameter. Hq is given in km/Mpc.s, 
A in units of c, to in Gy and Ui in units of oq. 

illustration of the importance of the marginalization process because it changes the peak values and credible 
regions depending if two or one-dimensional parameter space is used. The region of highest probability, at 
2cr level, is concentrated around a = with a large dispersion for A, as the curvature and/or the pressureless 
matter is fixed, this region displaces itself to positive values of a and to higher values of A, the opposite that 
happens with the PDF value of A, due in the last case to the marginalization procedure. The dispersion in 
A remains always large. It is obvious, on the other hand, from the two-dimensional credible regions that the 
case a = 1 is not ruled out. Using a larger sample of 253 data, the authors of reference |2H1 have constrained 
the parameters a and A such that [1.6,3.6] and [0.856,0.946] at 68% confidence level. But they have used 
the statistics and quartessence, and in this sense their result may be considered as compatible with ours. 
In reference ^51, restricting to the flat case, and fixing a = 1, for example, the authors have found that 
0.66 < A < 1 at 95% confidence level. This result is compatible with ours, see table 5. 

3.3 Estimation of Q^o 

The most general case where all five parameters are free in the GCGM predict fimo — O.OOIJ qo: while for 
the CGM, Qrno = 0.00+ 

o.oo- Hence, the unified scenario, with no pressureless matter, is favoured, but the 
dispersion is very high. Among the particular situations, where some of the parameters are fixed, the only 
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Table 5: The estimated parameters for the traditional Chaplygin gas model (CGM) and some specific cases 
of spatial section and matter content. We use the Bayesian analysis to obtain the peak of the one-dimcnsional 
marginal probability and the 2cr credible region for each parameter. Hq is given in km/ Mpc.s, A in units of 
c, ^0 ill Gy and in units of oq. 

relevant case is the flat Universe, which favours still the unified model, but with a much smaller dispersion. 
Repeating the analysis for ACDM, we find Vl^o — I-Olij gs. This prediction is a very important distinctive 
feature of ACDM with respect to the GCGM. In reference using a large sample of supernovae, the 
estimated value for the pressureless matter component, for ACDM, was fimo = 0.28 ±0.05 for a flat Universe. 
A similar value has been obtained in reference |17j . who employed also a very large sample of supernovae: 
0.29^0 03. In this case, our analysis gives fimo = 0.309lg'Q72, with a good agreement. The analysis of the 
WMAP data leads to Vlmo — 0.14 ± 0.02 -4 . AH these results are in some sense compatible due to the large 
error bar in our estimation. 

3.4 Estimation of Vl^q 

Remarkably, the GCGM, the CGM and the ACDM leads to almost the same predictions concerning the dark 
energy component when all parameters are free: 1.36to'85, l-34^o;7o and 1.36lo^7|, respectively. The best 
value does change appreciably when one or two parameters are fixed. The most important distinguishing 
feature is the dispersion, that is considerably higher, mainly in the upper uncertainty, for the GCGM. In 
comparison with the studies performed with a restricted sample of supernovae jl5lll6| . the main modification 
is the narrowing of the dispersion. The joint probability for Jlco and fi^o reflects what has been said in this 
and in the preceding sub-section. For the ACDM case, the two-dimensional picture for ri^o and ri„iO display 
an ellipse, with the best value around (1, 1), while the ellipse is highly distorted for the GCGM and CGM 
cases, with the best value around (1, 0). Moreover, the Icr, 2cr and 3cr contours are expressively displaced to 
the top of the diagram, mainly for the GCGM, what is a consequence of the high dispersion in this case. 
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Table 6: The estimated parameters for the ACDM model and some specific cases of spatial section and matter 
content. We use the Bayesian analysis to obtain the peak of the one-dimensional marginal probability and 
the 2(7 credible region for each parameter. Hq is given in km/Mpc.s, to in Gy and in units of ag. 

3.5 Estimation of Q^o 

Both in the GCGM and in the CGM, a closed Universe is clearly favoured: Qko = ^0-77^5 g4 and Qko ~ 
— 2.73^0 07, respectively. The dispersion, however is quite large. The quantity of pressureless matter displaces 
the maximum of probability to a value near or equal to zero and narrows the dispersion. Curiously, the 
dispersion for the GCGM has increased in comparison with the restricted sample of supernovae used in 
reference jJB], while for the CGM it has remained almost the same. A similar behaviour occurs in the 
dispersion of the other mass density parameters. For the GCGM the probability to have a closed Universe is 
96.29%; when the pressureless matter is fixed equal zero or equal to the baryonic content, this probability is 
smaller, 69.43% and 74.40%, respectively. The probability to have a flat Universe is 37.04%. For the CGM, 
such numbers increase slightly : for example, the probability to have a closed Universe is extremely near 
100%: 3.59(7. In this case, the probabihty to have a flat Universe is only 0.04%, but increases to near 60% 
after setting the pressureless matter. The ACDM model favours also a closed Universe: flko = ~2.12:ti'6i' 
but an open Universe is favoured when the pressureless component is flxed. Moreover, the probability to 
have a closed Universe in the ACDM is 98.08% when all free parameters are taken into account, but this 
value drops to zero when the curvature or the pressureless matter component is fixed. In principle the ACDM 
gives estimations closer to the CGM, but it must be remarked that the dispersion is very large. 

3.6 Estimation of Hq 

The predicted value of the Hubble constant today Hq is the most robust one. This can be understood by 
looking at the expression for the luminosity distance: the Hubble constant appears as an overall multiplicative 
factor. The predictions are: Hq = 65.00+};^^ for the GCGM and = 65.0lti;?i for the CGM. Fixing the 
curvature and/or the pressureless matter changes very slightly these predictions. Note that the dispersion 
is relatively small. In comparison with the restricted sample of references ^1 the best value for the 
Hubble constant has increased a little, and the dispersion has diminished. Repeating the analysis for the 
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Figure 2: The graphics of the joint PDF as function of (a, A) for the generaUzed Chaplygin gas model. The joint PDF peak 
is shown by the large dot, the credible regions of 1 cr (68, 27%) by the red dotted line, the 2 cr (95, 45%) in blue dashed line and 
the 3 (T (99, 73%) in green dashed-dotted line. The cases for QmO = are not shown here because they are similar to the ones 
with OmO = 0.04. 



ADCM, we find Hq = 65.00ti:74, very near the values found for the GCGM and CGM. The estimated value 
differs, on the other hand, from the predictions coming from CMBR, without a superposition of error bars: 
the WMAP predicts Ho^72± 0.5 j4j. 

3.7 Estimation of the age of the Universe, to 

The predicted age of the Universe for the GCGM is to = UMtllt Gy and for the CGM to = 12.73l:i:^4 Gy. 
These values are dangerously near the recent estimations age of the globular clusters to = 12.6^2 4 Gy. 
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Figure 3: The PDF of a for the generaUzed Chaplygin gas model. The solid lines are the PDF, the la (68.27%) regions are 
delimited by red dotted lines and the 2a (95.45%) credible regions are given by blue dashed lines. The cases for Qmo = are 
not shown here because they are similar to the ones with QrnO = 0.04. 



However, the error bars remain quite large. The estimation of the age of the Universe using the WMAP 
data gives 13.4 ± 0.3 Gy 4 . Compared with the previous analysis with a sub-sample of the supernovae, 
the predicted age has considerably diminished jl5lll(i| . Fixing the curvature and/or the pressureless matter 
increases slightly the predicted age (this is the opposite behaviour when sub-sample of 26 supernovae is 
used). Considering now the ACDM, the estimation of the age of the Universe leads to to = 12.70l^'22, ™ 
good agreement with the predictions of the GCGM and CGM. In the ACDM, however, values as high as 
17 Gy can be found by, for example, fixing ilmo — 0.04. Note that, for the GCGM, the best value for the 
product Hoto is 0.84, which is essentially the same for the other cases, while in reference 18 , for the ACDM 
model, Hoto = 0.96. 



3.8 Estimation of the deceleration parameter qo 

The value for the deceleration parameter qo in the GCGM with five free parameters is given by go = 
— 0.818lg'45g. The particular cases where the curvature and/or pressureless matter are fixed change very 
slightly this value. On the other hand, fixing a = 1, qo = — 0.883lg'429- These best values have not changed 
appreciably with respect to the restricted sample of 26 supernovae. However, the dispersion is considerably 
smaller. Repeating the analysis for the ACDM, we find qo — — 0.864lg;4Q5, barely differing from the previous 
models. In all cases, the probability to have a negative value for the deceleration parameter is equal to or 
very near 100%. 
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Figure 4: The one-dimensional PDF of A for the generaUzed Chaplygin gas model. The solid lines are the PDF, the \a 
(68.27%) regions are delimited by red dotted lines, the 2a (95.45%) credible regions are given by blue dashed lines and the 
3cr (99.73%) regions are delimited by green dashed-dottcd lines. The cases for Omo = are not shown here because they are 
similar to the ones with OmO = 0.04. 



3.9 Estimation of the moment the Universe begins to accelerate 

Another useful quantity is the redshift at which the Universe begins to accelerate, z;. Due to computational 
reasons, it is more practical to evaluate the value of the scale factor at the moment the Universe begins to 
accelerate a^, keeping in mind that the scale factor is normalized with its present value equal to unity. For 
the GCGM with five parameters, we find = 0.746;1;q'5i5- Imposing the curvature and/or the pressureless 
matter component does not change appreciably this result. For the CGM we find a, = 0.739l!5 o88- In 
reference the same analysis has been made for the GCGM with the restricted sample of supernovae: the 
value of a.i was smaller, but the dispersion was considerably higher. For the ACDM we have Oi = 0.732l°;?^i^. 
In all cases, the probability the Universe begins to accelerate before today is essentially 100%. All these 
results must be compared with that obtained in reference [l? ! which gives, translated in our notation, 
a. = 0.68t°;°^. 



4 Conclusions 

The aim of the present work was to present the most general analysis of the GCGM and CGM in what 
concerns the comparison of theoretical predictions with the type la supernovae data, using the 157 data of 
the "gold sample". All free parameters for each model were considered. In the case of the GCGM there 
are five free parameter: the Hubble constant Rq\ the equation of state parameter a\ the "sound velocity" 
related parameter A\ the curvature density of the Universe rjfco; the density parameter for the pressureless 
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Figure 5: The one-dimensional PDF of A for the traditional Chaplygin gas model (CGM). The solid lines arc the PDF, the 
la (68.27%) regions are delimited by red dotted lines, the 2a (95.45%) credible regions are given by blue dashed lines and the 
3(T (99.73%) regions are delimited by green dashed-dotted lines. The cases for OmO = are not shown here because they are 
similar to the ones with OmO = 0.04. 



matter (alternatively, the Chaplygin gas) fimo (alternatively, fico)- For the CGM, the number of parameters 
reduce to four, since a = 1. We have considered also the ACDM, where the number of parameters reduce 
to three: Hq, fifco and fimo (alternatively, fico)- 

A Bayesian statistical analysis was employed in order to obtain the predictions, with the error bars, for 
each parameter in each model. In doing so, the marginalization procedure was used so the one parameter 
estimations become more robust. This procedure consists in integrating in the remaining parameters in order 
to obtain a prediction for a given parameter. The results are the following. For the GCGM, we obtained: 
a = -0.75t^;°t = 65.00^^:^^, = = O.OOtJ;^^, = IMtlf^, ^ = 1.0001°;™°. 

The results for the CGM are: Hq ^ 65.0ll};^J, Vtm = -2.73tjj?, = O.OOtJ;^, ^co = l-34l[!;?o, 
A = 1.000t°;°?°. For the ACDM, we found: Hq = &bmt\ll, ^ka = -2.12^^;^^, = I-Oltass' ^co = 
LSG^Q ^g. We have also evaluated the age of the Universe, the value of the decelerating parameter today 
and the moment the Universe begins to accelerated. The results for the GCGM, the CGM and the ACDM 
are, respectively: to = 12.63t? ?^, 12.7it\\\, 12.70t?-^^; qo = -O.SlStjjl^, -0.883l|]-^^g, -0.864t°-^°°; 
a. = 0.7461H?!, 0.739t°:°^«, 0.732:H?i 

In the GCGM, the best value for a is negative but, due to the large dispersion, high positive values are 
also allowed. This may be compared with the results of references [2Hll2niEnii were a takes high positive 
values. In reference 23 j the apparent discrepancy is due to the quartessence choice and the statistical 
method employed : the authors used the statistics in order to obtain the confidence regions. Other 
crucial parameter is A. Both for the GCGM and CGM, the best value is in principle 1, but the finite 
resolution used in the numerical computation suggests that the peak value of A can be between 0.98 (or 
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Fi gure 6: The graphics of the joint PDF as function of (Q^o, ^co) for the generalized Chaplygin gas model, where p{QmOy f^co) 
is a integral of p{a, Hq, f2mOi ^cOi ^) over the (a, Hq,A) parameter space. It shows a different shape for the confidence regions 
with respect to the KCDM but similar to the CGM of ref. 1161 . The joint PDF peak has the value 0.584 for (n^Oif^co) = 
(0.00,1.00) (shown by the large dot), the credible regions of 1 cr (68,27%, shown in red dotted line), 2cr (95,45%, in blue 
dashed line) and 3(t (99,73%, in green dashed-dotted line) have PDF levels of 0.104, 0.037 and 0.010, respectively. As 
f^feO + ^mO + f^cO = li the probability for a spatially flat Universe is on the line Q,mO + f2cO = li above it we have the region 
for a closed Universe (fc > 0, H^q < 0), and below, the region for an open Universe (k < 0, Ofeo > 0)- 

0.99) and 1.00. Futhcrmore, the particular cases of fixed curvature and matter densities shows a value near 
but small than 1 as the best value of A, such that the ACDM case [A — 1) is almost ruled out. 

The results reported above indicate that, for the GCGM, CGM and ACDM, a closed Universe is favoured. 
For the GCGM and CGM the unified scenario (quartessence) , where the pressureless matter density is 
essentially zero, is also favoured. In any case, a small fraction of pressureless matter must be introduced 
in order to take into account the baryons. However, the dispersions are quite large, decreasing significantly 
for the flat Universe case. It is curious that in the ACDM, the density parameter for pressureless matter is 
around unity. The predictions for the dark energy component are similar for the three models, but in the 
GCGM the dispersion is very expressive, allowing for very large values for the density parameter of the dark 
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Fi gure 1 '. The graphics of the joint PDF as function of (ilmOi^co) for the traditional Chaplygin gas model (CGM), where 
p(£7mOif^co) is a integral of p(-ffoi ^mOi f^cOi ^) over the {Ht},A) parameter space. It shows a different shape for the confidence 
regions with respect to the ACDM but similar to the CGM of ref. 1161 . The joint PDF peak has the value 1.473 for 
(HmO, flco) = (0.00, 1.22) (shown by the large dot), the credible regions of 1 cr (68, 27%, shown in red dotted line), 2 u (95, 45%, 
in blue dashed line) and 3ct (99, 73%, in green dashed-dotted line) have PDF levels of 0.542, 0.103 and 0.009, respectively. As 
f^feO + ^mO + f^cO = li the probability for a spatially flat Universe is on the line HmO + f^cO = li above it we have the region 
for a closed Universe (fc > 0, < 0), and below, the region for an open Universe (k < 0, f2feo > 0)- 

energy. 

The analyses made predict a not very old Universe for all three cases. The best value is compatible, 
for example, with the estimated age of the globular clusters However, it must be stressed that such 

compatibility is verified mainly because of the high dispersion in those evaluations. The predictions for the 
deceleration parameter and for the moment the Universe begins to accelerate does not vary expressively in 
the diS'erent models studied. 

In general, the predictions above agree with those obtained by refs. ^|^], where the supernovae data 
have been studied extensively in the context of the ACDM. In what concerns the previous studies of the 
GCGM using the supernovae Cni~EZIi there are many differences due mainly to the range of values assumed 
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Figure 8: The graphics of the joint PDF as function of (f2mO,f^co) for the ACDM model model, where p(nmOif^co) is a 
integral of p{Ho, QmO, ^co) over the Hq parameter space. The joint PDF peak has the value 4.22 for (QrnO, ^co) = (1-01, 1.35) 
(shown by the large dot), the credible regions of 1 (68, 27%, shown in red dotted line), 2 a (95, 45%, in blue dashed line) and 
3cr (99, 73%, in green dashed-dotted Une) have PDF levels of 1.380, 0.247 and 0.025, respectively. As Uko + ^mO + ^cO = 1, 
the probability for a spatially flat Universe is on the line QmO + f^cO = 1, above it we have the region for a closed Universe 
{k > 0, QicO < 0)i below, the region for an open Universe (k < 0, QkO > 0). 

for the parameter a: the authors considered < a < 1, the statistical methods are different, and frequently 
some parameters were considered fixed. In general, those previous results obtained a positive value with a 
great degeneracy for a ,24 . However, by allowing a to take large positive and negative values, the conclusions 
of those works may somewhat change. Ref. |30| is in agreement, as the the GCGM and CGM are preferred 
over the ACDM, and the flat Universe case is also favoured (like ref. [211 suggests) in most GCGM and 
CGM cases. 

In references [151 a more selective sample of 26 supernovae has been used. In comparison with the 
present results, the main differences are the following: the predicted value for a has become slightly less 
negative; the predictions for A remained essentially the same, except for the CGM, where it became nearer 
unity; the value of the Hubble parameter has increased slightly; the age of the Universe has become smaller; 
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Figure 9: The one-dimensional PDF of QfcOi f^mO a^nd f^cO for the generahzed Chaplygin gas model. The solid lines are the 
PDF, the la (68.27%) regions are delimited by red dotted lines, the 2cr (95.45%) credible regions are given by blue dashed 
lines and the Scr (99.73%) regions are delimited by green dashed-dotted lines. As QcO = 1 — f2feo ~ ^feOi for ^mO = we have 
^cO = 1 — f^feOi foi' f^mO = 0.04 then QcO = 0.96 — QkO '^'^'i for Q^o = we also have QcO = 1 — f^mO- 



the value of the deceleration parameter today has changed slightly. In general, the dispersion has (sometimes 
only marginally) decreased with the large sample of supernovae, except for the parameters fifeo, f^mo smd 
flco where the dispersions were, for the GCGM, almost always smaller with the restricted sample of 26 SNe 
la. This behaviour may be an issue for future analyses of the GCGM and CGM using thousands of SNe la 
(from SNAP and other projects). 

The traditional Chaplygin gas model, where a = 1, remains in general competitive. When the five 
parameters are considered, the probability to have this value is of 19.58%. But, this probability increases as 
much as to 72.70% if the curvature is fixed to zero and if only baryons account to the pressureless matter. 
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Figure 10: The one-dimensional PDF of QkOt ^mO 3.nd QcO for the traditional Chaplygin gas model. The solid lines are the 
PDF, the 1(T (68.27%) regions are delimited by red dotted lines, the 2(T (95.45%) credible regions are given by blue dashed 
lines and the 3cr (99.73%) regions are delimited by green dashed-dotted lines. As QcO = 1 — f^feo ~ f^feOi for f2mO = we have 
f^cO = 1 — f^feOi for OmO = 0.04 then QcO = 0.96 — QkO s-nd for Q^o = we also have QcO = 1 — ^mO- 

So, as far as the type la supernovae data are considered, the Chaplygin gas scenario is not ruled out. 

As a logical future step, we plan to cross the estimations from different observational data, like gravita- 
tional lensing, the large scale structures data (2dFRGS), the anisotropy of the Cosmic Microwave Background 
Radiation (CMBR) and X-ray gas mass fraction data EHl EH EHl ESI ■ 
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PDF of Qi,o for A-1 pop of Q^o for A=l 
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Fi gure 111 The one-dimensional PDF of ^ko^ ^mO and for the ACDM model. The solid lines are the PDF, the Icr 
(68.27%) regions are delimited by red dotted lines, the 2a (95.45%) credible regions are given by blue dashed lines and the 3(t 
(99.73%) regions are delimited by green dashed-dotted lines. As ilcO = 1 — f^fco ~ f^feo, for ^mO = we have QcO = 1 — f^fto, for 
S^mO = 0.04 then QcO = 0.96 — QkO E'nd for Q^o = we also have QcO = 1 ^ ^mO- 
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